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This research article presents the non-linear dynamic of a three-link robotic 
manipulator formulated by the Newton-Euler method. The planar 
manipulator is composed of three links and three revolute joints rotating 
about the z-axis. The three nonlinear non-homogeneous dynamic equations 
have been solved graphically with the help of MATLAB by phase variable 


method. The work represents the graphical solution of the transient response 


of angular position, and angular velocity of each link member for a 
Keywords: predetermined interval of time. With the help of simulated value from 
MATLAB, torque characteristics have been determined for different torque 


Genetic algorithm ratios and optimum torque has been derived using a genetic algorithm to 


MATLAB move the manipulator in a proper direction. 
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1. INTRODUCTION 

Robot manipulators can be assumed to be a chain of link mechanisms with highly nonlinear 
dynamics. The forward and inverse kinematics are presented using the D-H convention and the deduced 
mathematical model has been inferred by Newton-Euler approach. To obtain the optimal performance of the 
robot manipulator, a precise dynamic model of the robot manipulator is required. Garg ef al. [1] determined an 
optimal path using a genetic algorithm and optimization was achieved using simulated annealing (SA). 
Exhaustive simulation was conducted for different types of manipulators namely MELFA RV-14A having six 
degrees of freedom and three planar revolute joints to find the distance between the end effector and the object 
using an artificial neural network [2]. The payload of the mobile robotic arm is determined by a nonlinear 
control law using optimal feedback [3] designed for a given trajectory task. This law is given by the solution 
by the iterative method of a sequence of nonlinear Hamilton-Jacobi-Bellman equations. Aghanouri et al. [4] 
developed a manipulator actuated by DC motors. The optimum path is derived by optimizing system 
parameters to minimize performance indices including energy [5]. Talezadeh et al. [6] presented the nonlinear 
dynamics of two-link manipulators for dynamic modeling of the system by the Lagrange equation of motion 
where optimum control was adopted to analyse the motion of the manipulator. Some literature has presented 
the whole-body dynamics of nonlinear equations of the hybrid cable-driven robots (HCDRs) where new 
methods were developed to solve redundancy. Compared to the existing methods, torque optimization for 
actuated and un-actuated joints can solve resolution problems as well as active satisfactory disturbance [7]. 

Baressi et al. [8] describe the kinematic modeling using the Denavit Hartenberg (DH) convention 
and the dynamic modeling robotic arm has been presented by Lagrange-Euler (LE) and Newton-Euler 
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algorithms. Different algorithms such as artificial intelligence, genetic algorithm, simulated annealing, and 
differential evolution have been adopted and compared to provide the best results for minimization of torque. 
Televnoy et al. [9] provides Lagrange equations of motion of a nonlinear matrix equation. The experiment 
has been performed to find out the solution of a moving object between two points. The angular displacement, 
speed, and acceleration for six links of the manipulator have been presented. Agustian et al. [10] present the 
vision-based robotic manipulator where inverse kinematics using pseudo-inverse Jacobian (PIJ) and DH 
forward kinematics have been adopted and controlled by a proportional derivative controller. The sorting task 
depending on color was made to evaluate the error to implement the manipulator on a real system. The reviews 
was given [11] on dynamic analysis and intelligent control techniques for flexible robot manipulators. A 
comparative study of dynamic analysis and control strategies was presented for flexible manipulators. 
Korayem et al. [12] have determined the non-linear dynamics and control of flexible mobile manipulators 
focusing on the determination of maximum payload. Wu et al. [13] present the minimum actuator torque 
range by torque optimization of a 3-DOF parallel manipulator. Two approaches, Lagrangian and primal-dual 
neural network, were presented together to make real-time optimization of joint torque for kinematically 
redundant manipulators [14], [15]. Gao et al. [16] proposed joint torque optimization of flexible manipulators 
with redundant space and vibration suppression where the Lagrange method has been adopted to represent 
the dynamics of the robot arm. Naghshineh and Keshmiri [17] investigated an over actuated system for 
dynamic cost function by applying real time optimization on a cooperative robot. Wolniakowski et al. [18] 
presented a method required for joint torque minimization of a serial manipulator to determine the optimal 
task placement. Woolfrey et al. [19] described a control strategy for minimizing joint torque using null space 
control of a redundant manipulator where the dynamic torque has been reduced by applying an external force 
to the gripper element. Agbaraji et al. [20] presented the design of the manipulator by calculating and 
analyzing the joint torques by evaluating performance in terms of speed and displacement of the arm based 
on the predetermined values of the torques. Singh et al. [21] find the position vectors of the six-arm robot by 
forward kinematics and joint angles by inverse kinematics in MATLAB with the help of a robotic toolbox. 
Brandstotter et al. [22] and Petrenko et al. [23] presented a generalized closed-form solution of the dynamic 
models of parallel robots using some simple Jacobian matrices where the dynamics of the legs were expressed 
in the joint coordinates and that of the platform in the form of cartesian variables. Sun et al. [24] proposed the 
methodology for deriving the closed-form inverse kinematic solutions of the 6-DOF robot on the position 
level where the analytical inverse solution of all the joints was given out and compared with that of the forward 
kinematic solution. The mathematical formulation and simulation of the two-link planar robot manipulator 
with forward kinematics and dynamics were presented with the help of the D-H convention and Newton-Euler 
method where the results in terms of joint angles and angular velocities have been presented graphically [25]. 
In this paper, the Newton-Euler formulation has been adopted to present the dynamics of the three- 
link manipulator. The solution relating torque and angular displacement is obtained analytically by expressing 
the three second-order non-linear non-homogeneous differential equations into six first-order differential 
equations by phase variable method. The angular positions and speed of each link have been presented 
graphically for constant torques of fixed duration using MATLAB. Then optimization is done using a genetic 
algorithm to obtain minimum torque that can be applied to each link-joint actuator to obtain optimum results. 


2. THREELINK ROBOT MANIPULATOR 

The length of the three links of the robotic manipulator is given by l, L; and l5. Three joints are 
named J,,Jz,J3, respectively, as in Figure 1. mi, m», ms are the mass of the first link, second link and third 
link, respectively. Link parameters are considered to be a; = a; = «4 = 0. All the three joints are revolute 
joints rotating about the z-axis, and all the links are rigid links. The angular rotations of each link-joint 
combination are denoted as 01,02, 03. 

The initial conditions are assumed to be 


Wy 20,09 =0,% =0 


where Wo, initial angular speed, «y , initial angular acceleration, Vo, initial linear velocity and Vg, initial linear 
acceleration. Therefore, 


ù = [0 9.81 oO]? 
qi= [61.05.05]; 

qi= 191, 82, 83] 

Gi = (01, 02, 03]; 

Link variable: F;, fi, nj, tj 
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where F; is the total external force exerted at the center of i^ link, f; is the force exerted on i link by i-1* 
link, n; is the moment acting on i™ link by i-1* link, 7; is the torque on i" joint. 


Figure 1. Coordinate assignment of the robot manipulator 


3. NEWTON-EULER FORMULATION FOR COMPUTATION OF JOINT TORQUE 

Computational work has been carried out using Newton — Euler methodology for transformations of 
three link coordinates to formulate the speeds of three-link manipulator. Initial speed of each link is assumed 
to be zero. For the kinematic model, computation of each joint and link variables is required for computation 
of each joint torque to move the manipulator in a desired direction. Forward kinematics and backward 
kinematics of motion has been applied in the Newton — Euler approach as presented in (1) to (3). T1,T2,T3are 
the joint torque applied to the joint actuator for link i = 1, 2, 3, respectively. 


Tz = imjl,l, [cos( 6; + 05), + sin( 0, + 05); ^] + malls [sin 6; (6, + 62)? + 
cos 63 (6, + 6)] + = Mls" (B, + 6, + 63) + zmslg cos(0, + 0; + 04) (1) 


Tz = 5Mgl,l,[cos( 2 + O3)6, + sin(6; + 03)0?] + : msl,ls [sin 85 (8, + 82)? + 

cos 04 (6, + 6)] + im BË, + 6, + 63) + =M3l3g cos(@, + 0; + 63) + (Mm; + 
zmjll;[cos 6; Ë, + sin 03 02] + (m; + m;)l2 (8, + 62) + msl,l [cos 0; (0, + 

6, + Öz) — sin 0; (6, + 62 + 63)"] + (mz + 2 m;)L,g cos( 6, + 02) (2) 


7, = 2msl, l5 [cos( 0 + 058, + sin( 02 + 63)6, ^] 2 malls [sin 05 (6, + 62)” + 

cos 03 (6, + 07)] + Em $Ë, + Ö, + Öz) + ~m3l3g cos( 6, + 05 + 04) + (m3 + 

zm3)ll, [cos 6; Ë, + sin 8; 02] + (ms  2m;)l2 (B, + 82) + 2 mslola[cos 05 (Ö, + 

6, + 63) — sin 64 (0, + 0; + 0,7] + (m; + ~m)log cos( 0, + 05) + (m3 + m; + 

=m,)126, + (mz + mj)lil; [cos 6; (6, + Öz) — sin 6; (6, + 6)?] — 

<mslyls[sin( 02 + 05) (0, + Öz + 63)? — cos( 02 + 03) (6; + Öz + Ö3)] + (ms + 

m, + =m lig cos 0, (3) 


4. COMPUTED JOINT TORQUE CHARACTERISTICS OF EACH LINK-JOINT PAIR 

The dynamic equations of motion of each link-joint coupling are highly non-linear in nature as shown 
in (1) to (3). To solve the above equations of motion relating t;and 0; (for i" link), phase variable model has 
been adopted. The three non-linear non-homogeneous 2™ order differential equations are converted into six 
first-order differential equations by phase variable method. Assuming Z4, Z2, Z3, Z4, Z5 and Z, are six phase 
variables where: 


24 = 6, (4) 
z = 6, (5) 
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Z4 = 8, (6) 
Z4 = 01 (7) 
Zs = 6, (8) 
Ze = 0, (9) 
Therefore: 
4-8 (10) 
ż = Ö, (11) 
255p. (12) 
ż4 = EE (13) 
ig =O, =z (14) 
Ż6 = 0,27, (15) 


Non-dimensionalizing (1), (2) and (3), we have 


T3 — ~cos(@, + 63)6, + zsin(6, + 0,)6,* + <sin 05 (0, + Å>)? + Z cos 05 (6, + 


6.) += (6, + Öz + Öz) + =cos( 0; + 05 + 03) (16) 


ie = cos( 6; + 63)6, + =sin( 0, + 04)82 + zsin 63 (0, + 62)? + = cos 63 (6, + 
6,) += (6, + Öz + Öz)  2cos (0, + 0z + 04) + =cos 0; 0, - sin 0, 02 += (6, + 
6.) + ~cos 63 (6, + 0; + öz) — zsin 03 (0, + 0; + 63)? + = cos( 6, + 8) (17) 


Rm ~cos( 6; + 04)8, + zsin(6; + 03), + zsin 63 (0, + 05)? + = cos 63 (6, + 
65) + - (6; + 6, + 63) + ~cos( O, + 0, + 03) + = cos 0, 0, + sin 0,02 + = (6, + 

65) + ~cos 6, (0, + 8, + 63) — <sin 0, (0, + Ô, + 63)? + = cos( O, +0) + “6, + 

= cos 6, (0, + 2) — = sin 6, (0, + 04)? — zsin(6, + 04,)(0, + Åz +63)? + 

~cos( 0, + 03)(6; + Öz + Öz) += cos 6, (18) 


Using (16) to (18), 6,, 6 and Ë; can be expressed as a function of 0,, 05, 93 and their derivatives. Therefore, 
six first order equations can be formed from (10) to (15) and can be solved using MATLAB which gives the 
relationship of angular displacement and applied torque for various torque ratios as shown in Figures 2 to 8. 


Angular Rotation Vs. Applied Torque 
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Figure 2. Angular position versus applied torque (11: 12: 13: 1:1:1 (in Nm); time: 0.5sec) 
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Angular Rotation Vs. Applied Torque 
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Figure 3. Angular position versus applied torque (11: 12: 13: 2.25:1.5:1 (in Nm); time: 0.5sec) 
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Figure 4. Angular position versus applied torque (11: 12: 13: 4:2:1 (in Nm); time: 0.5sec) 


Angular Rotation Vs. Applied Torque 
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Figure 5. Angular position versus applied torque (11: 12: 13: 6.25:2.5:1 (in Nm); time: 0.5sec) 


Angular Rotation Vs. Applied Torque 
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Figure 6. Angular position versus applied torque (11: 12: 13: 9:3:1 (in Nm); time: 0.5sec) 
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Angular Rotation (rad) 


Angular Rotation Vs. Applied Torque 
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Figure 7. Angular position versus applied torque (11: 12: 13: 16:4:1 (in Nm); time: 0.5sec) 
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Figure 8. Angular position versus applied torque (11: 12: 13: 25:5:1 (in Nm); time: 0.5sec) 


5. TORQUE OPTIMIZATION USING GENETIC ALGORITHM IN MATLAB 

To optimize 13, joint torque in the third joint actuator has been plotted against t2 and 05 as shown in 
Figure 9. Genetic algorithm has been applied to minimize the objective function given in (19). Optimized 
value of 04 has been evaluated by plotting 05 with respect to evaluated 04and 05 as shown in Figure 10. The 
optimized results of all the required parameters are given in Table 1. 


Theta2 


Figure 9. Surface generation of T3 for torque 


* 3 vs. T2, Theta2 


i 0.5 
en £ T2 Theta2 AD 30 Thetat 


optimization and 02 


Table 1. Optimized torque and corresponding angular displacement 
ti(in Nm). tin Nm) _13(inNm)  0;(rad) (rad) _03(rad) 
31.25 12.5 1.3658 1.1376  -0.1861 1.3202 


165 


è  Theta3 vs. Theta, Thet: 


Figure 10. Surface generation of 05 with respect to 0; 
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74 = a + bsin(mnx,x;) + ce- r2» (19) 


where a, b, c, m, and w are constant coefficients of the objective function (19) and evaluated in MATLAB as 


a — 2.683 
b — 0.9224 
c — —0.4006 
m = 0.2138 
w = —0.6463 


x,and xz are the torque required and resulting angular displacement of the link2, respectively. Boundary 
condition has been placed as 1 € x4 € 12.5 and —3 € x; <0. 
From Table 1, it is evident that the results show a successful optimization using genetic algorithm on 


the problem of torque optimization. For practical implementation, assumptions can be made that 0 < 04 < = 
0z0,- 7 and 0 < 04 < 7 and the condition|6,| < 20,must be followed if 0; is positive and 05 found to be 


negative. Therefore, the result shown is to satisfy the above-mentioned condition and the joint torques as well 
as corresponding angular displacement may be evaluated as given in Table 1. 


6. CONCLUSION 

The manipulator dynamics, as stated in (2) to (4), represent the non-linearity of the robot arm 
system. Here the open chain of links and its corresponding behavior of the manipulator has been studied 
which shows the angular position vs. applied torque characteristics. The optimum value of torque with 
respect to the desired movement of the manipulator arm has been evaluated using a genetic algorithm with 
boundary conditions in MATLAB which also satisfies the relative angular position constraints. Also, the 
system may be assumed to be a series of three inverted pendulums exhibiting the nature of a chaotic system 
which is certainly nonlinear in nature. Therefore, stability of the system may be achieved with the help of a 
sliding mode controller as well as by behavior-based control where it is to split a complex dynamic into 
several simple equations which are quietly related to the problem stated in this research paper. 
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